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Abstract 

We present a very promising method for a measurement of the Higgs boson parity 
using the H/A — > t + t~ — » p + v T p~v T — > it + it v t it~ it v T decay chain. The method 
is both model independent and independent of the Higgs production mechanism. 
Angular distributions of the r decay products which are sensitive to the Higgs 
boson parity are defined and are found to be measurable using typical properties of 
a future detector for an e + e~ linear collider. The prospects for the measurement 
of the parity of a Higgs boson with a mass of 120 GeV are quantified for the case 
of e + e~ collisons of 500 GeV center of mass energy with an integrated luminosity 
of 500/& -1 . The Standard Model Higgsstrahlung production process is used as an 
example. 
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1 Introduction 



If the Higgs mechanism is realized in nature, a complete determination of its properties is 
among the central tasks for future colliders. In particular, determining the CP properties 
of the Higgs boson is one of the important goals of a future e + e~ linear collider operating 
at a center-of-mass energy between 350—1000 GeV 0-0] • Different models, such as the 
Standard Model (SM) , the general Two Higgs Doublet Models (2HVM), the Minimal 
Supersymmetric Standard Model (AiSSAi), and others predict different properties for 
the details of the Higgs boson signature. It is of great importance to verify that, whatever 
scenario is realized in nature, the scientific program of a linear collider will be able to 
distinguish among these models. 

Several methods have been proposed to distinguish between a scalar (J vc = ++ ) and 
a pseudoscalar (J vc = _+ ) spin zero Higgs particle PHlO"! ■ In this study we investigate 



the case of a Higgs boson which is light enough that the W + W~ decay channel remains 
closed. Then the most promising decay channel for the model independent parity deter- 
mination is H/A — ► t + t~ . It was previously proposed M that the best r decay channel 
for parity determination would be — > tc^u. The weakness of this method is that the 
Higgs boson rest frame needs to be precisely reconstructed. 

To address resolution issues it is necessary to perform Monte Carlo studies where 
the significant details of theoretical effects and detector conditions can be included. In 



Ref. [1.1 J we have extended the algorithm of Refs. [12, 13] for generating Higgs decay 
(independently from its production mechanism) to the pair of r-leptons including the 
complete spin correlation matrix used in the subsequent decay of the r leptons. The 
reaction chain e + e~ — > Z(H/A), H/A — > t + t~; r ± — > ■k ± u t {u t ) was studied. It was found 
that even small effects of smearing seriously deteriorate the measurement resolution. An 



independent study [14|, using the PANDORA Monte Carlo generator of Ref. [H| confirms 



that result. This leaves no doubt that if we want to use the Higgs boson r decay channel 
for the measurement of the Higgs parity, a better technique is required. 

In this paper we continue to investigate Higgs decay into a t + t~ pair. We extend our 
study to the H — > t + t~ decay with r 1 * 1 — ► p ± v T (v T ) (i.e. the channel with the largest 
branching ratio) and then p ± — > tt^tt . Spin correlations are distributed over three levels 
of the decay chain. Complicated geometrical distributions need to be defined. We have 
found the Monte Carlo method particularly useful already at the level of defining the 
observables. All of the main results are cross-checked using two independent analyses 
based on an extension of universal interface (Ref. [fLl}| ) for TAUOLA, and with the PANDORA 
Monte Carlo generator |0| which is interfaced to PYTHIA 6. 1 [ll| for fragmentation and 
decay processes. 

The rest of the paper is organized as follows. The theoretical considerations which are 
used to understand the decay chain of the Higgs particle are explained in section 2 where 
we also give some details of the TAUOLA Monte Carlo simulation. In section 3 we define 
the observable we use to distinguish between the scalar and pseudoscalar Higgs boson. 
In section 4 we list our assumptions on smearing and we discuss imposed cuts, detector 
effects and necessary adaptations introduced into our observables. Main numerical results 
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are also given in this section. Section 5, a summary, closes the paper. 



2 Determination of the CV quantum numbers of the 
Higgs boson 

The H/A parity information must be extracted from the correlations between r + and r _ 
spin components which are further reflected in correlations between the r decay products 
in the plane transverse to the r + r~ axes. This is because the decay probability, see |4J], 

T{H/A^r + r-) ~ 1 - s\ + s\~ ± s T ± s T ± ~ (1) 

is sensitive to the r^ 1 polarization vectors s T and s T+ (defined in their respective rest 
frames). The symbols ||/_L denote components parallel/transverse to the Higgs boson 
momentum as seen from the respective rest frames. This suggests that the experi- 
mentally clean t + t~ final state may be the proper instrument to study the parity of the 
Higgs boson. 

The spin of the r lepton is not directly observable but it manifests itself in the distri- 
butions of its decay products. Depending on the decay channel, the polarimetric strength 
is different. The first Monte Carlo program for e + e~ colliders where the density matrix of 



the r lepton pair was used was KORALB [17,18]. Let us recall some basic properties of that 



solution, which we have adopted to the case of the Higgs boson decay in Ref. [JTTJj . The 
algorithm is organized in two steps. In the first step, the r lepton pair is generated and 
the r leptons are decayed in their respective rest frames as if there were no spin effects 
at all. In the second step, the spin weight is calculated and rejection is performed. If the 
event is rejected, only the generation of the r lepton decays is repeated. The spin weight 
is given by the following formula 



1(1+ E ''V''/'') (2) 



wt 
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where, as a consequence of formula ([TJ), the components i? 3j3 = —1, R^i = ±1, i? 2 ,2 = ±1 
(respectively for scalar and pseudoscalar) and all other components are zeroQ. 

In the following, we focus on the — > p ± u decay channel. It is interesting because 
it has, by far, the largest branching ratio (25%). However, in comparison to r ± — » ix^v 
decay, its polarimetric force is more than a factor of 2 smaller. It was found, see e.g. [ 19 , 20 1 , 
that in many cases this can be improved if information on the p decay products, i.e. on 
details of the decay r ± — > 7r =t 7r°zy, are used. This is of no surprise because the polarimetric 
vector is given by the formula 

h l = Af(2(q-N)q l -q 2 N l ^j (3) 



1 See Rcf. |yj for detailed definition of the quantization frames used for the spins of the t + and r 
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where M is a normalization function, q is the difference of the 7r ± and 7r° four-momenta 
and N is the four-momentum of the r neutrino (all defined in the r rest frame) see, 
e.g. [21]. Obviously, any control on the vector q can be advantageous. It is of interest to 
note that in the r lepton rest frame, when m n ± = m n o is assumed, the term 

q-N= (E n ± - E n o)m T . (4) 

Thus, to exploit this part of the polarimetric vector, we need to have some handle on the 
difference of the ir^ and 7r° energies in their respective r leptons rest frames. Otherwise, 
the effect of this part of the polarimetric vector cancels out and one is left with the part 
proportional to the p (equivalently v) momentum. Without using the energy difference 
we would arrive at a case nearly identical to the one where the r decays to a single ir 
except with the disadvantage that the coefficient diminishes the spin effects by more than 
a factor of 2. 

Already from this preliminary discussion we realize that the appropriate observable 
must rely on a distribution constructed out of at least 4 four-vectors. The Monte Carlo 
method is already essential at the level of designing the observable. Our search for im- 
provements of the results obtained in began with a study of acollinearity distri- 



butions defined in the Higgs boson rest frame, when instead of the r — > nv decay mode 
the mode r — ► pv was used. The difference between the scalar and pseudoscalar Higgs 
particle which was still visible for r — > ttv turned out to be practically invisible in the 
case of r — > pv once detector smearings were introduced. We found these results rather 
discouraging, and we have turned our attention to another possibility. As inspiration we 
have used two methods, one used at LEP 1 for the measurement of the r polarization and 



another one, proposed for low energy e + e colliders, see e.g. |L9 



2.1 The Monte Carlo 

For the following discussion all the Monte Carlo samples have been generated with the 
TAUOLA library 0-0]. The PHOTOS f2|,[25| Monte Carlo program could be used for 



generating radiative corrections in the decays of the Higgs boson and r leptons, but it 
was switched off. For the production of the r lepton pairs the Monte Carlo program 
PYTHIA 6. 1 |16|] is used.0 The production process e + e~ — > ZH — > p + pT (qq)H has been 



chosen with a Higgs boson mass of 120 GeV and a center-of-mass energy of 500 GeV. 
The effects of initial state bremsstrahlung were included in the PYTHIA generation. 

For the r lepton pair decay with full spin effects included in the H — > t + t~, r ± — > 
p ± v T (v T ), p ± — > 7T ± 7r° chain, the interface explained in Ref. [O was used. It is an extended 
version of the standard universal interface of Refs. []i"2~, 13]|. For the sake of confidence we 



have confirmed all numerical results presented in this paper with the second simulation 
using the PANDORA Monte Carlo generator fllS] . 



2 It was shown that the interface can work as well in the same manner with the HERWIG [Em generator. 
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3 The acoplanarity of the p + and p decay products 

In this section we advocate a new observable where we ignore the part of the polarimetric 
vector proportional to the p (equivalently v) momentum in the r rest frame. We rely 
only on the part of the vector due to the differences of the ^ and ir° momenta, which 
manifests the spin state of the p ± . 

In the Higgs rest frame the p momentum represents a larger fraction of the Higgs's 
energy than the neutrino. Therefore, we abandon the reconstruction of the Higgs rest 
frame and instead we use the p + p~ rest frame which has the advantage that it is built 
only from directly visible decay products of the p + and .0 

We take for both p's the n^ir decay channel. In the rest frame of the p + p~ system 
we define the acoplanarity angle, (p*, between the two planes spanned by the immediate 
decay products (the 7r ± and 7r°) of the two p's. 

3.1 Defining an optimal variable 

The variable <p* alone does not distinguish the scalar and psuedoscalar Higgs. To do 
this we must go further. The r — > tt ± tt°i' spin sensitivity is proportional to the energy 
difference of the charged and neutral pion (in the r rest frame), see formula (|4]). We have 
to separate events into two zones, C and V, 



C : 2/iy 2 > 

V : y x y 2 < 



where, 



^ _ E w + — E w o _ ^ _ E n - — E T o 
E^+ + E n o E v - + E^o 



E T ± and E w o are the vr* 1 , 7r° energies in the respective r ± rest frames. If the cuts are applied 
using the four-momenta from the generation level boosted to the r rest frame without 
any smearing then we call them respectively Chare and T>b ar e- If the cuts are applied using 
the smeared four- momenta boosted to the replacement r ± rest frames (defined below), 
then they are called C reco and V reco . 

In Fig. [I] we plot the distribution of (p*, where the left hand plot contains the events 
where the energy difference between the 7r + and 7r° defined in the r + rest frame is of 

3 The use of correlation angles to measure the Higgs parity have already been proposed, see, e.g. (2(]] 
or ||, but their definitions depended on the r and Higgs (or Z/j) rest frames which are difficult to 
measure. In our approach, these frames are replaced by others which are easily measured and yet retain 
significant sensitivity to the Higgs parity. In some special cases our approach may open the way to 
studying the Higgs parity at hadron machines such as the Tevatron or the LHC. 



4 



17.1 ■ 10- 3 
8.6 ■ 10- 3 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

Lf Lf 

Figure 1: The p + p~ decay products 7 acoplanarity distribution angle, (p*, in the rest frame 
of the p + p~ pair. A cut on the differences of the tt^ tt° energies defined in their respective 
rest frames to be of the same sign, selection Cb are , is used in the left plot and the 
opposite sign, selection T>^ are , is used for the right plot. No smearing is done. Thick 
lines denote the case of the scalar Higgs boson and thin lines the pseudoscalar one. Units 
valid for the 500 GeV e + e~ CMS (scalar 120 GeV mass) Higgsstrahlung production only. 
Otherwise arbitrary units. 

the same sign as the energy difference of tc~ and tc° defined in r~ rest frame (selection 
Chare and formula ffy. The right hand plot contains the events with the opposite signs for 
the two energy differences (selection Vb are ). It can be seen that the differences between 
the scalar and pseudoscalar Higgs boson are large. If the energy difference cut was not 
applied, we would have completely lost sensitivity to the Higgs boson parity. 

Unfortunately, since the r-lepton is not measurable, such a selection cut cannot be 
used directly. We now go on to define our choice for detection parameters. Then we will 
define realizable methods for making the energy cuts, C reco and V reco , and we will discuss 
phenomenologically sound results. 

4 Detector Effects 

To test the feasibility of the measurement, some assumptions about the detector effects 
have to be made. We include, as the most critical for our discussion, effects due to 
inaccuracies in the measurements of the ir momenta. We assume Gaussian spreads of 
the 'measured' quantities with respect to the generated ones, and we use the following 
algorithm to reconstruct the energies of 7r's in a measurable substitute for their respective 



rest frames. 



1. Charged pion momentum: We assume a 0.1% spread on its energy and direction. 



achieved with a 15% energy error and a 27r/1800 direction error in the gammas 



3. The replacement r lepton rest frames and the n and ir energy differ- 
ences: To make the measurement we replace the difficult to measure Higgs rest 
frame with the p + p~ rest frame. Define replacement four-monenta, p±, in the p + p~ 
rest frame for the unmeasurable four-momenta: 

a) Define p° + = p°_ = mu/2. 

b) Define the direction of the p + and p_ three-momenta to be the direction of the 
p + and p_ three-momenta in the p + p~ rest frame, respectively. 

c) Define the magnitude of the p± three-momenta so that p\ = p 2 _ = m 2 . 

Boost the 7r + , 7T°, 7r~, 7T° momenta to the respective rest frames, p + and p_, of their 
replacement r + and r _ . The 7r energies defined this way are used in the C reco and 
V reco energy difference cuts. 

When we compare predictions for scalar and pseudoscalar, we should consider not 
only properties of its decay, but we should take care of the possible differences in the 
production mechanisms as well. To avoid multitude of options we have excluded this 
point from our study. We use Higgsstrahlung production mechanism both for the scalar 
and pseudoscalar Higgs boson. Except for the size of the cross section, our analysis and 
its conclusions do not depend on the choice. We have checked that it is indeed the case 
by varying beam energies and choosing other production mechanisms Q. 

4.1 The results of detector smearing 

If we use the true, generator level, r rest frame to define the energy differences between 
the p decay products but smear the momenta of the pions when used in the calculation of 
the acoplanarity angle, (p°, the resulting distributions (not shown) are very similar to the 
unsmeared case of Fig. [|. When we use the selection cuts C reco and V reco (and thereby use 

4 We have studied a W/Si sampling EMCal with 1.25 m inner radius and projective towers with 1800 
segments in both 9 and 4> in a 5 Tesla field. The direction of reconstructed photons can be defined by 
the center of energy of the hit cells in a cluster defined by contiguous hit cells. The combination of 
the small Moliere radius of tungsten, the fine segmentation and using detailed hit information yields a 
direction resolution of about 1/6 27r/1800 which is well below the resolution assumed above, see e.g. |27f| 
and references therein for details. 

5 Note that large Higgsstrahlung cross-section arise only if Higgs boson has a sizable scalar component. 
In such a case our method could measure its pseudoscalar admixture. 
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the replacement r rest frames as well as smearing the 7r momenta) we obtain the results 
shown in Fig. 0. We see that the effects to be measured diminish but remain clearly 
visible []. 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 0.0 0.5 1.0 1.5 2.0 2.5 3.0 



cp' cp' 

Figure 2: The p + p~ decay products' acoplanarity distribution angle, if', in the rest frame 
of the p + p~ pair. A cut on the differences of the tt^ tt° energies defined in their respective 
replacement r ± rest frames to be of the same sign, se is used in the left plot 

and the opposite sign, selection T> reco , is used for the right plot. All smearing is included. 
Thick lines denote the case of the scalar Higgs boson and thin lines the pseudoscalar one. 
Units valid for the 500 GeV e + e~ CMS (scalar 120 GeV mass) Higgsstrahlung production 
only. Otherwise arbitrary units. 



4.2 The potential measurement resolution 

To determine the Higgs parity at an operating linear collider, a set of events will be col- 
lected over a period of time. An event selection will be made from the data set to isolate 
Higgsstrahlung events where the process H/A — > t + t~ — > p + v T p~v T — > 7t + 7t v t 7t~7t°i> t 
occurred. This sample will be reconstructed and the distribution of the measured variable, 
ip', defined in sections 3 and 5, will be compared with simulated reconstructed distribu- 
tions for a scalar and a pseudoscalar Higgs such as those shown in Fig. ||]. A goodness of 
fit test such as an unbinned maximum likelihood will be performed for both hypotheses. 

6 We have studied several options for the definition of the separation cuts C reco and D reco . In one case 
we have directly used the smeared laboratory frame energies. In another, we have used all the information 
available from the reconstruction of the Higgs boson rest frame. All of these choices lead to practically 
identical versions of Fig. @. 
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One hypothesis will be favored over the other and statistical techniques will be applied 
to estimate the confidence level of the conclusion. 

For present purposes, it is important to note that the level of certainty obtained will 
depend upon the specific sample collected. Another run under the same conditions with 
the same integrated luminosity would, in general, result in a sample of a different size and 
with different fit results due to the random sampling effects inherent in all measurements. 
Thus, we cannot a priori predict the level of certainty that will be achieved when the 
measurement is eventually made. However, we can simulate a large set of possible data 
samples that might be collected and gather some insight into the range of certainties 
possible. 

We have applied a binned maximum likelihood technique to a set of approximately 
1300 simulated data samples where half the samples were derived from a CV even Higgs 
and half from a CV odd Higgs. The Higgs mass was taken to be 120 GeV, the beam energy 
was 250 GeV per beam. The p^ were reconstructed from their gammas and ir ± decay 
products. The 7r ± were smeared 0.1% in energy and direction, the gamma energy was 
smeared 15% and their direction by 27r/1800, a typical calorimeter cell size. We assumed 
an integrated luminosity of 500 fb" 1 . Beamsstrahlung and ISR effects are included. To 
account for detector acceptance effects, event selection efficiency and impurity we assumed 
an overall efficiency of 60%. Studies |28|,|29| have shown these are realistic estimates of 
event selection efficiencies and purities. 

Based on these assumptions, using a binned likelihood fit and comparing each data 
sample with the distributions for a CV even and a CV odd Higgs and comparing the 
resulting fits, we find that every data sample will identify the correct parent distribution 
with a confidence level of at least 95% and 86% of all samples will make the correct 
identification above the 3a confidence level. Thus, we see that we have an excellent 
chance of correctly determining the CV of the Higgs with this technique. 

The technique is quite robust relative to the measurement resolutions of the charged 
pions and the gammas (from the neutral pions). For example, decreasing the direction 
resolution of the gammas even by a factor up to six has a negligible effect on the overall 
measurement resolutions cited above. The (p* plots show that the odd and even states 
have on average large angular differences, thus small errors in particle resolutions will not 
change if* of an event enough, to make a significant impact on the overall distribution. 



5 Summary 

We have studied the possibility of distinguishing a scalar from a pseudoscalar couplings of 
light Higgs to fermions using its decay to a pair of r leptons and their subsequent decays 
to r ± — > u T p ± and p^ — > 7r ± 7r°. 

We have discussed an observable which is very promising. Using reasonable assump- 
tions about the SM. production cross section and about the measurement resolutions 
we find that with 500 of luminosity at a 500 GeV e + e~ linear collider the CV of 
a 120 GeV Higgs can be measured to a confidence level greater than 95%. To confirm 
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the method we have used two distinct Monte Carlo programs and |TT1,|T3H and the 
observables were coded independently. We emphasize that the technique is both model 
independent and independent of the Higgs production mechanism and depends only on 
good measurements of the Higgs decay products. Thus, this method may be applicable to 
other production modes including those available at proton colliders as well as at electron 
colliders. 

Finally we note that several improvements are still possible. If, instead of unweighted 
events in our distributions as given in Fig. ^|we use events weighted with weight = [ 2/x2/a | 
(see formula [5]) the statistical significance will increase^ by a factor of about 1.5. Other 
multi-meson final states such as r — > 7t + 7t + tt~i> can be used to increase the statistical sam- 
ples. As argued in |Kj they should lead to additional data of spin significance comparable 
to our t — > pv channel (if the appropriate observables are found). Finally, we may expect 
that a measurement of the r flight direction may turn out to be helpful, similarly as it was 
the case for the r polarization measurement at LEP 1 |3l],|3^]. For each new production 
mode that can be analyzed using this technique the sample size increases accordingly. 
We may expect the final statistical significance for the parity measurement to be at least 
factor of 2-3 times better than in our conclusions. Such improvements will lead to more 
complex observables with many special cases and require good control of the systematic 
errors. We leave them to future studies but we think that, because of such possibilities, 
our estimation of the parity resolution will turn out to be largely conservative. 
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